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THEORETICAL AND EXPERIMENTAL INVESTIGATION OF HEAT 
CONDUCTION IN AIR, INCLUDING EFFECTS OF OXYGEN DISSOCIATION 

By O. I red k rick Hansen, Hichard A. Early, Frederick. K. Alzofon, and Fred C. Witte born 


SUMMARY 

The one-dimensional, non steady flow of heat 
through gases at constant pressure is considered 
where the coefficients of thermal conductivity and 
di flu si city are functions of temperature . For ideal 
gases, the di if u s i city is taken to he directly propor- 
tional to the integral of thermal conductivity. In 
the case of air pit is further assumed that the diffusiv- 
ity becomes a constant in the temperature interval 
where oxygen dissociation occurs. These assump- 
tions facilitate numerical solution of the heat flow 
equation for the case of a semi-infinite gas medium 
with a constant initial temperature and a constant 
bo undary te mpe rature . 7 'he solutions are used , 
along with measured shock velocity and the measured 
heat flux through a surface from, which the shock wave 
reflects , in order to evaluate the integral of thermal 
conductivity as a function of temperature. This 
integral replaces temperature as the dependent vari- 
able for heat conduction in media which have a vari- 
able coefficient of thermal conductivity. The experi- 
mental results agree with theoretical estimates 
within the accuracy of the measurements arid within 
the uncertainty in theory . 

INTRODUCTION 

In the classical solutions for heat transfer by 
conduction, it has usually been assumed that the 
coefficients of thermal conductivity and of diffusiv- 
ity are constant (see refs. 1 and 2, for example). 
This assumption is quite good for many solids and 
liquids if the absolute temperature throughout the 
medium does not vary by large factors. In gases, 
on the other hand, the coefficient of thermal con- 
ductivity is approximately proportional to the 
half power of temperat ure if the gas is inert (ref. 3). 
For constant pressure processes, the diffusivity of 


an inert gas is approximately proportional to the 
three-halves power of temperature. Even greater 
variations with temperature occur in gases which 
are chemically reactive. Thus, the classical heat- 
conduction solutions apply to gases only where the 
differences in absolute temperature throughout the 
gas are very small. This condition is not fulfilled 
in many practical heat-conduction problems, of 
course. For example, it has been shown that, 
neglecting radiation, the heat transfer occurs es- 
sentially by a conduction process at the stagnation 
region of vehicles traveling at high speed through 
t ho atmosphere (ref. 4). Within this same region, 
the temperature, and therefore the other prop- 
erties of the gas, changes by large factors. This 
heat-conduction problem is further complicated 
by the fact that the air molecules may dissociate in 
the flow around vehicles flying faster than about 
7,000 feet per second. In addition, appreciable 
ionization may occur at speeds above 30,000 feet 
per second. Such chemical reactions absorb heat 
in regions where dissociation and ionization occur 
and liberate heat again in regions of recombina- 
tion, thus behaving like heat sinks and heat sources 
in the gas, respectively. Ncrnst (ref. 5) shows 
that such effects can strongly influence the heat- 
conduction characteristics of the gas. 

The steady-state heat-conduction problem can 
be characterized in terms of the temperature and 
the coefficient of thermal conductivity, provided 
reaction rates are rapid enough to insure chemical 
equilibrium locally throughout the gas. Hirsch- 
f elder (ref. 6) has analyzed the one-dimensional, 
steady-state heat conduction in a chemically re- 
acting gas in which large variations in the coeffi- 
cient of thermal conductivity appear. Coffin and 
O'Neal have experimentally verified Hirsclifelder’s 
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results for the ease of dissociating N 2 0 4 (ref. 7). 
In the ease of nonsteadv heat conduction the co- 
efficient of thermal difFusivity also appears. It is 
the purpose of this paper to develop solutions for 
one-dimensional, nonsteady heat conduction in 
gases at constant pressure that take into account 
variations in both the coefficients of thermal con- 
ductivity and diffusivitv. These solutions will he 
used in conjunction with experimental measure- 
ments in a shock tube for the purpose of evaluating 
an average coefficient of thermal conductivity as 
a function of temperature. 
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SYMBOLS 

coefficient of thermal diffusivitv, ~y f - 

P c p 

stoichiometric coefficient for component 
Ju in a chemical reaction (eq. (21)) 
reacting cluanical component 
specific heat per unit mass 
binary mass diffusion coefficient for 
diffusion of particles type l into 
particles type j 

coefficient of thermal conductivity 
coefficient k averaged over tempera- 



chemical equilibrium constant in terms 
of partial pressures 

exponent of temperature giving the 
variation of k with temperature, also 
the number of mols of gas per unit 
volume 
pressure 

heat flux per unit area 
gas constant per unit mass 
time 

absolute temperature of gas medium 
temperature at which the coefficient of 
thermal conductivity is a maximum 
temperature at the interface between 
gas and solid media, also used as a 
standard reference condition 
initial temperature of solid medium 

dimensionless variable. ?/-*/— n - 
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distance 

mol fraction of the it h component of a 
gas mixture 
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dimensionless variable, 

a 4«n t 

dimensionless integral of thermal con- 
duct iv it v, “ 

reference value for the parameter z 
(see eq. (lh)) 

reference for the limiting value of z 
(see eq. (H))) 

coefficient of resistance change with 
tempera tu re- 

density 

integral of thermal conductivity, 



value of the parameter <p foi* normal 



gradient operator 

2 p , 

— I e-^du, 

7T Jo 

SUBSCRIPTS 

conditions at temperature 7 T 0 , also the 
condition at the interface between 
gas and solid media 
initial conditions in the shock tube 
conditions following the initial shock 
wave produced in the shock tube 
conditions following the reflection of 
the initial shock wave 
initial conditions in the shock-tube 
reservoir 

initial conditions in the gas medium, 
identified with condition 3 iu the 
shock-tube experiments 
conditions at temperature T c 
conditions in the gas state at constant 
pressure 

conditions in tin* solid state medium 
bounding the gas phase 


THEORETICAL CONSIDERATIONS 

The problem which will be considered is the 
nonsteadv, one-dimensional heat conduction 
through a semi-infinite gas medium at constant 
pressure with uniform initial conditions. First, 
the analysis will treat the conduction of heat 
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through an ideal gas to a constant temperature 
reservoir. Next, the same problem will be con- 
sidered for a gas which is chemically reactive, 
such as dissociating air. Finally, the analysis 
will treat the heat conduction through the gas to 
a semi-infinite solid which also has uniform initial 
conditions. 

The partial differential equation which describes 
nonsteady heat flow by conduction is the familiar 
diffusion equation 

^>-^-V(KT)=° (1) 


pendent variable in these solutions. This func- 
tion is the product of the temperature and an 
average coefficient of thermal conductivity (eq. 
(2)). It may be noted that the lower limit of the 
averaging integral is arbitrary, just as the reference 
level normally used for temperature is arbitrary. 
Absolute zero is used as the lower limit here as a 
convenience. The heat llux at any position is 
just the gradient of the <p function at that point 
in the medium. 

2 = “~grad ^ (5) 

HEAT CONDUCTION THROUGH AN IDEAL CAS 


The temperature T is the dependent variable, C*, 
is the specific heat per unit mass, and p is the 
density of the medium. Usually, the coefficient, 
of thermal conductivity k is considered constant 
so that it can be taken outside the differential 
operator. However, for the problem of interest 
here, conduction of heat through gases, k must he 
treated as a variable. In such problems it is 
generally permissible to treat the coefficient k as a 
function of temperature only, however. Then a 
useful transformation is 

*>=( kdTsslcT (2) 

and the differential equation (1) can be put into 
the form 

-^-«VV= 0 (3) 


In the ease of gases, the differential equation (J) 
must he solved with a variable coefficient of 
diffusivity. For ideal gases the approximation is 
made that k varies as the 7/th power of tempera- 
ture and that the density p is given by the ideal 
gas law. Then t he diffusivity is 


lik o 


T n j 1 


( 6 ) 



ko 

1)7 


_7 T n + l 
0 


(7) 


where the subscript zero refers to values at some 
standard condition. Thus, diffusivity is propor- 
tional to <p 




( 8 ) 


where a is the thermal diffusivity 


It may he noted that equation (2) is the same- as 
the Kirchoff transformation referred to by Jakob 
(ref. 2), and has also been used by Langmuir (ref. 
8) in calculating the rate of heat loss from solid 
surfaces immersed in a gas. (-rank (ref. 9) uses 
a similar transformation in treating some mass 
diffusion problems. 

It is apparent from equation (J) that all of the 
heat-conduction solutions which have been worked 
out on the basis of a constant, thermal conductivity 
apply equally well to media for which thermal 
conductivity varies with temperature if the thermal 
diffusivity is relatively constant. The function <p 
merely takes the place of temperature as the de- 


This relation is reliable over a considerable range 
of temperature 4 for monatomic gases and is a rea- 
sonably good approximation for many polyatomic 
gases, including common diatomic gases, such as 
hydrogen, oxygen, and nitrogen. Sometimes a 
better approximation is given by a linear relation 
with a finite intercept, but this can always be' 
transformed to a direct proportionality by appro- 
priate adjustment of the lower limit for the integral 
(p. In the integrations which follow, it will be 
necessary to use only the result that the diffusivity 
a is approximately proportional to <p. 

Consider now the nonsteady, one-dimensional 
flow of heat into a constant temperature reservoir 
from a semi-infinite medium. The heat -conduc- 
tion equation reduces to 

^-a-^=0 (9) 

d 1 or 
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Boltzmann (see ref. 9 or 10) lias shown that for 
certain boundary conditions, and provided that 
the coefficient of diffusivity is a function of the 
dependent variable only, the solution to equation 
(9) may be expressed as a function of x f\t. This 
occurs when the boundary conditions in both 
time and space can also be expressed as functions 
of x/\t alone, so that x and t arc not involved 
separately. In the present case, the diffusivity is 
taken to be a function of <p alone, and the boundary 
conditions are constant, independent of x and t . 
Therefore the solution may be given as a function 
of the dimensionless time-distance parameter 


V 



( 10 ) 


It will also be convenient to make the dependent 
variable dimensionless with the transformation 


z=v!<Po ( 11 ) 

Then the partial differential equation (9) becomes 
the ordinary differential equation 


For convenience in plotting, the variable y is 
replaced by w 



so tint 1 :- z will approach its limit in each case at 
about the same point, that is where w is ap- 
proximately 2.0. It is also convenient to use the 
function. 

s*=l + (2„ — 1) erf w (16) 

as a normalizing parameter so that the solutions 
can be plotted accurately on a single graph for a 
wide range of conditions. The parameter z* is 

TABLE I. — BOUNDARY CONDITIONS IN A SEMI- 
INFINITE GAS MEDIUM WITH A CONSTANT 
INITIAL TEMPERATURE, r«>, AND A CON- 
STANT BOUNDARY TEMPERATURE, T 0 

. A '”' 

I k (IT 


a d~Z 

(in <hf 


+-’/ 



( 12 ) 


V 

\ 'W 


The factor 4 under the radical of equation (10) 
is an arbitrary stretching factor chosen so that the 
dependent variable reduces to its usual form in 
the case where diffusivity is constant. The stand- 
ard state will he taken as the constant conditions 
which exist at the boundary between the conduct- 
ing medium and the reservoir at £=(). Then the 
boundary conditions are 

*o=l (13a) 

Zoz=<p*!*pq (13b) 


'Where a is proportional to <p , 
becomes 


d 2 z. ky dz _ 
, >T ~i ~ =0 

dy dy 


equation 


( 12 ) 

(14) 


It does not seem promising to seek an exact 
analytic solution to equation (14) which will 
satisfy the given initial and boundary conditions. 
However, if one specifies the value of (dz/dy ) 0 , it 
becomes relatively easy to integrate the equation 
numerically. Some results of integrations of 
equation (14) are shown in figure 1 and table I. 
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just the function to which 3 reduces in the case 
where thermal diffusivity is constant. 

Figure 1(a) shows the ratio z/z* for specific 
values of the limit 3^, which have been obtained 
bv numerical integration of equation (14). The 
initial slopes corresponding to these values of z m 
may be found in table I. The ratio z/z* is every- 
where larger than unity and goes through a 
distinct maximum. The effect is more pro- 
nounced the larger the value of The integral 
<P, however, is a smooth, monotonically increasing 
function of w or ?/, similar to the error function in 
appearance, but rising more steeply near the 
origin as indicated by the ratios z/z*. The 
integral may be calculated from 

(17) 


From equation (5), the heat flux through the gas 
is 


< Po dz 
\ 4 a 0 t dy 


(18a) 


or in terms of the normalizing function z* 


___<P0 /dz \ 
\n(r x t \dz*) 



(18b) 


Figure 1(b) shows the ratio dzjdz* as a function of 
w for various values of z m . The curves show how 
the heat flux varies throughout a gas in comparison 
with the first derivative of the error function. 

The first derivative of z may be thought of as 
a normalized dimensionless heat flux (see eq. 
18(a)). The value of this normalized heat flux 
at the boundary determines the value of (see 
table I). The quantity z m * will be defined as 
the limiting value which z would approach, for the 
given boundary derivative, if the thermal diffus- 
ivity were constant; that is, 



Figure 1(c) shows the ratio 2 «>/£«,* as a function 
of (dz/dy) 0 . It will be noted that z m /z m * ap- 
proaches a linear function of (dz/dy) 0 as the 
derivative inc r eases . 

HEAT CONDUCTION THROUGH CHEMICALLY REACTIVE GASES 

The results obtained so far are for inert gases. 
Where the gas is chemically active and in local 


equilibrium, Hirschfeldcr (ref. G) shows that the 
heat transferred by diffusion of reacting molecules 
and atoms results in increased thermal conductiv- 
ity. If the lien t of reaction is large, the coefficient 
of thermal conductivity becomes much larger 
than normal and goes through a pronounced 
maximum in the interval of temperature where 
the chemical reaction occurs. The thermal dif- 
fusivity, on the other hand, is not greatly affected 
by the reaction because the increase in thermal 
conductivity is compensated hy a proportional 
increase in the specific heat As a result, the 
integral <p increases very rapidly with temperature 
compared to the diffusivity. Considered as a 
function of <p f the diffusivity is approximately 
constant during the interval of temperature where 
the coefficient of thermal conductivity is greatly 
increased due to the chemical reaction. 

The diffusivity is most nearly independent of <p 
when* increases most rapidly with temperature, 
that is, where the coefficient of thermal conductiv- 
ity is a maximum. The temperature T c where 
this maximum occurs is approximately the 
temperature where the specific heat is a maximum, 
or where the equilibrium mol fractions of the 
reacting species change most rapidly with temper- 
ature. The value of T c may be calculated more 
precisely from equations derived by Butler and 
Brokaw for t he coefficient of thermo 1 conductivity 
(ref. 11) 




w — -~ 

ji nx ( D Sj 


(b i.i j bjj /) 


( 20 ) 


where K p is the chemical equilibrium constant in 
terms of partial pressures, n is the number of mols 
per unit volume, the x t are the mol fractions of the 
ith component, and the 6* are the stoichiometric 
coefficients for the reacting components /?, when 
the chemical equation is written in the form 

S 6iB t ==0 (21) 

i 

For most common gases in the equilibrium stab', 
these quantities are all precisely known functions 
of temperature. The binary mass diffusion co- 
efficients Dij arc not so well known, in general, 
but they are relatively slowly varying functions 
of temperature and do not strongly influence 
the position of the maximum in k. Any reason- 
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(a) Integral of thermal conductivity as a fun *tion of distance and time. 


l'i<! r hk 1. One-dimensional heat flow through an ideal gas where ’he initial temperature and boundary temperature 

are constants. 
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(b Heat flux as a function of distance and time. 
Figure 1 . — Continued. 
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Dimensionless heat flux at the be undary, I 


(c) Limit for the integral of thermal conductivity as a function of the boundary heat flux. 

Figure 1. Concluded 


able approximation for the mass diffusion co- 
efficients will suffice for the calculation of T c . 

In the approximations which follow, it is 
assumed that a is proportional to <f> up to the 
temperature 7 T C , above which it is a constant a c . 
Once tlie boundary temperature 7 T 0 is chosen, the 
value of z at the transition point is given by 




The numerical integration of equation (14) is 
terminated at the point where z reaches z c . This 
point uniquely determines values for the argument 
y c and t fie derivative ( dz/dy ) c . Since the dif- 
fusivitv is assumed to he constant for z^>z c , the 
solution can he continued analytically to yield 


N Zc/\ \ z c \Z C / 


The value of z^ is easily found from equation (23) 
by let til g y go to infinity. Figure 2 shows the 
ratio z m f z m * for chemically reacting gases as a 
function of the derivative ( dz/dy ) a . Again, 
is given by equation (19). Curves are shown for 
equal to 10, 20, 30, 40, 50, and 60. One can 
interpolate between these curves for the specific 
value of z c given by equation (22), depending on 
the bom dary condition and the chemical reaction 
consideivd. The dashed curve shows the inert- 
gas solution from figure 1(c) for comparison. 

To illustrate an application of the above analy- 
sis, cons der the case of air in which oxygen dis- 
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Figcre 2. Integral of thermal conductivity its a function of boundary heat flux for a chemically reacting gas with large 

boat of reaction. 


sociation can occur. The temperature T c for this 
case is shown as a function of pressure in figure 3. 
The exponent n approaches a value of about 1/2 
for high temperature air. Then if the boundary 
temperature is taken to he 293° K, z c takes the 
values 21.3, 24.9, 29.6, 35.3, 42.2, 5 0.5, and 60.0 
for air pressures of 0.0001, 0.001, 0.01, 0.1, 1.0, 10, 
and 100 atmospheres, respectively. Oxygen dis- 
sociation occurs at lower temperatures when pres- 
sure is decreased. Thus, the lower the pressure, 
the smaller the value of z Cy and the smaller the 
value of which results iti a given heat flux at 
the boundary (see fig. 2). Similarly, is de- 
creased when the boundary temperature is in- 
creased (eq. (22)) and again smaller values of 
are required to produce a given slope (dzjdy) o. 

The functional relation between difFusivity and 
the integral <p is shown for air in figure 4. The 
ordinate and abscissa are normalized by the factors 
z c a 0 and z c <p () , respectively, so that the relation 
becomes independent of pressure. The solid line 
on the figure shows the function assumed in the 


integration of equation (12). The dashed curve 
shows a theoretical estimate of the actual function 
based on the transport coefficients for air which 
are given in reference 12. This relation does not 
depend strongly on numerical accuracy of the 
coefficients. Tlius, the exact function is probably 
not greatly different from the estimate, even 
though the latter is based on elementary kinetic 
theory and very approximate collision cross sec- 
tions. Jt will he noted that the theoretical dif- 
fusivity is somewhat smaller than the assumed 
value in the interval where dissociation first begins, 
hut is larger in the interval where oxygen dis- 
sociation nears completion. These errors tend to 
be averaged out by the integration process. 
When the oxygen dissociation becomes complete, 
the difFusivity again increases rapidly with <£>, and 
the present solution loses validity. The splitting 
of the dashed curve in figure 4 indicates that the 
normalization factors z c clq and z c <Po do not remove 
the pressure dependence of this point. 
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coefficients of thermal conductivity and diffusion 
and uniform initial conditions. Mersman, Berg- 
gren, and Boelter (ref. 13) show rigorously that 
the temperature at the interface between the two 
media comes instantly to a constant value, inde- 
pendem of time. It can be shown analytically that 
the saint 1 result applies where thermal conduc- 
tivity is an arbitrary function of temperature, pro- 
vided the thermal diffusivity is constant in each 
medium . In the case of gases, the diffusivity is not 
constant and the solutions are not in analytic form 
amenable to a general solution for the interface tem- 
perature. However, it can be seen that, a particu- 
lar constant value of the temperature 7 7 0 will just 
balance the heat flux between the two media for 
all values of time t. Note that the heat flux at 
the boundary of a gas, with constant boundary con- 
dition v'o? varies inversely as the square root of 
time, 


/<V\ = <Ai /dz\ 

\ d-iy o A 4a 0 1 \<ty ) o 


(24) 


Figure 3. — Pressure and temperature conditions in air at 
which oxygen dissociation is responsible for a maximum 
in the coefficient of thermal conductivity. 

HEAT CONDUCTION BETWEEN TWO SEMI-INFINITE MEDIA 
WITH UNIFORM INITIAL CONDITIONS 

The foregoing results apply to the heat conduc- 
tion through a semi-infinite gas medium with a 
constant boundary temperature. Such a condi- 
tion is approximately satisfied where the gas is 
bounded by a thick wall of a metal which has a 
very large thermal conductivity and heat capacity 
per unit volume compared to the gas. In this 
case, the temperature change in the metal is small 
compared to the temperature change through the 
thermal diffusion layer of the gas, and the metal 
behaves essentially as an infinite heat reservoir at 
nearly constant temperature. Thus, the solutions 
veil in the previous sections apply directly to 
such a case. 

It is important to inquire now whether the solu- 
tions for the constant boundary condition will 
also apply to the noiisteady heat conduction 
through a gas bounded by a medium with arbitrary 
characteristics. The second medium might be 
an insulating type material or another gas at a 
different initial temperature, for example. If is 
instructive to recall first the case of the sudden 
thermal contact between two media with constant 


since the derivative (< dzjdy ) 0 is a function only of 
(Pn and (see fig. 2), and is independent of time. 
By defii ition, the constant value of v? 0 implies a 
constant boundary temperature, 7 7 0 . Conse- 
quently, T 0 can always he adjusted to a particular 
constant value which will satisfy given uniform 
initial conditions in two adjoining semi-infinite 
media a; well as continuity of heat flux beween 
the two media. 

Consider now the ease where the wall is a solid 
having constant coefficients of thermal conduc- 
tivity k s and of diffusivity Equating the heat 
fluxes at the boundary, one obtains 



Normalized infegral of thermal conductivity, — 7 — 

z c 9 0 

1 igure 4. — Diffusivity for air as a function of the thermal 
conductivity integral. 
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<Po ( dz \ JcAT* T x ) ^ 25 ) 

y'4 a t) t \<ty / o \ 7r a v £ 

where jT s is the initial temperature of the solid. 
This result suggests a method of combining the 
analysis of nonsteady heat conduction through a 
gas with some experimental measurements for the 
purpose of determining the conductivity integral 
<p as a function of temperature. A shock wave 
reflected from the closed end of a tube can bo used 
to create a layer of gas with essentially zero How 
and constant, high-temperature initial conditions. 
Such an arrangement lias been used by Smiley, 
for example, for making measurements of thermal 
conductivity in high-temperature argon (ref. 14). 
Assuming equilibrium, one can calculate the tem- 
perature of the shock heated gas from measured 
shock velocity or pressure. Then, if the tempera- 
ture at the surface of the shock-tube wall is meas- 
ured, and the thermal coefficients of the wall 
material and of the gas at low temperature are 
known, the normalized heat flux at the boundary 
{(Izjdy) o may be evaluated from equation (25). 
This derivative uniquely determines the value of 
upon integration of equation (12), or alter- 
natively may be evaluated graphically from the 


solutions plotted in figure 2. Thus the integral 
is determined as a function of T m , the initial 
temperature of the shock heated gas. The subse- 
quent sections of this report will discuss such 
experiments performed to determine some values 
of <p for air. 

EXPERIMENT 

SHOCK-TUBE APPARATUS 

The shock tube which was used to produce 
liigh-temperature air for the heat-conduction 
measurements is illustrated in figure 5. The 
tube consists of a cylindrical steel pipe, 2 ^-inches 
inside diameter and 22 feet long. The reservoir 
chamber is separated from the shock tube by a 
copper diaphragm. The diaphragm is ruptured 
when desired by a spring driven punch. The 
pyramidal tip of the punch cuts the diaphragm 
so that it tears into four equal leaves which fold 
back against the tube as shown in figure 5(e). In 
order to assure a clean rupture it is necessary to 
choose the thickness of the diaphragm so that it 
is stressed to about 85 percent of its yield strength 
under the pressure of the reservoir gas. After 
the rupture, the reservoir gas flows into the tube 
forming a shock wave in the gas with which the 
tube is initially filled. 



(a) Photograph of .shock tube. 
Fioure 5 . - Shock-tube apparatus. 
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Helium is used for the* reservoir gas since it lias 
a low molecular weight and the maximum ratio of 
specific heats. Because of these properties, strong 
shock waves are produced with moderate reservoir 
pressures (see, for example, the shock-tube equa- 
tions in ref. 15). In addition, the reservoir is 
made large enough to function essentially as an 
infinite chamber in order to maximize the shock 
strengths obtained. This also serves to reduce 
the strength of those expansion waves, which 
propagate back into the reservoir, to a negligible 
amount. Thus, these expansion waves do not 
limit the interval of steady test conditions. 

After the shock wave reflects from the end of 
the tube, the high-temperature air produced is in 
a steady-state condition, except, for possible non- 
equilibrium effects. At the temperature and 
pressure conditions produced for the present 
experiments, the chemical relaxation times appear 
to be shorter than the available test interval (see 
appendix A), and it will be assumed here that the 
air is in equilibrium. The steady-state interval 
persists until disturbances arrive, several hundred 
microseconds later, due to t Ik* interaction between 
the reflected shock wave and the air-helium 
interface. 

A distance-time diagram of the wave patterns 
in the shock tube is shown in figure 5(c) for con- 
venience in establishing notation. The initial 
conditions are designated by a subscript 1 ; condi- 
tions following the initial shock wave arc 1 identified 
by subscript 2; and conditions following the 
reflection of the shock wave* from the end of the 
tube are denoted by subscript 3. The initial con- 
ditions in the reservoir are assigned the subscript 4. 


It is possible to calculate all of the equilibrium 
thermodynamic properties of the gases in regions 
2 and 3 as a function of the initial gas temperature 
and pressure, T x and p u and the speed ?/, with 
which the shock wave propagates into region 1. 
Table II lists the nominal values for 7\ and p lf 
which were used for the present tests, and the 
shock-wave Mach number Af l} which is obtained 
when helium at 54 atmospheres and normal tem- 
perature is used in the reservoir. The calculated 
pressures, temperatures, and densities of the air 
in region 5 (p 3 , 7 7 3 , and p 3 ) are given. These are 
taken from the graphs presented in reference* l(j. 
The table also shows the calculated interval of 
steady conditions in region 5 and the relaxation 
times for tin* vibrational excitation of nitrogen 
molecules and for the dissociation of oxygen 
molecules. These arc discussed further in 
appendix A. 

CALIBRATION MEASUREMENTS 

Tin* initial pressure of air in the shock tube, p h 
was measured with a resolution of about iO.] 
mm Hg. The initial temperature, T t , was 
measured within 0.5° 0. The shock Mach 
number, Af u was determined by measuring the 
time interval required for the shock wave to 
travel between two pressure detectors mounted 
18 incl es apart in the side wall near the end of 
the slack tube (stations 1 and 5 shown in figs. 
5(b) a id 5(d)). The pressure detectors were 
piezoelectric barium tilanalc crystals embedded 
in an e >oxy resin. The electrical output, from the 
crystals was used to start and stop an electronic 
chronograph having a resolution of 0.1 miero- 


TARLJ: II. -NOMINAL SHOCK TURK OPERATING CONDITIONS: INITIAL TEMPERATURE, T x - 295° K; 
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(c) 



(d) 

(b Schematic diagram of shock tube. 

(c) Time-distance wave diagram for shock tube flow. 

(d) Schematic drawing of shock tube test section. 


Figure 5. — Continued. 
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the side wall mounted gage are shown for a num- 
ber of different shock speeds in figures 0(a) through 
0(c). The pressure jumps from jh to p 2 when the 
initial shock wave passes, and then later jumps 
to ]h when the reflected shock arrives. (See 
sketch (a).) 



At higher shock strengths, some disturbances 
arrive at the side wall position ahead of the 
reflected shock, apparently through the boundary 
layer. This effect can be seen in figure 6(c), for 
example. In order to avoid this effect above 
Mach number 5.5, the pressure p z was measured 
with the detector mounted at station 0 in the end 
cap of the shock tube. A larger quartz piezo- 
electric pressure detector with a 50 ke resonant 
frequency was used at this station. The spurious 
side wall effects were absent at station 6 as shown 
in figures 6(d) and 6(e). The oscillations on all 
the pressure records art 1 due to excitation of the 
lowest resonant frequency of the pressure pickups. 
These have been averaged out in determining the 
pressures from the records. Due to this factor 
and the other irregularities, uncertainties in 
measured pressures are the order of 10 percent. 

Figure 7(a) shows the pressure p 2 as a function 
of the theoretical p> corresponding to the measured 
Mach number M v . The agreement is reasonably 
satisfactory in view of the uncertainty in the data. 
Figure 7(b) shows the measured pressure ^ 3 as a 
function of the calculated values for p 3 . Again 
the agreement is within the order of accuracy of 
the measurements. The value of /> 3 is consider- 
ably more sensitive to nonequilibrium effects than 
is p 2 (see ref. 16), and so the comparison shown in 
figure 7(b) indicates that the air reaches a state 
close to equilibrium. 


TEMPERATURE MEASUREMENTS 

A temperature detecting instrument was 
mounted flush with the wall in the end cap of the 
shock tube (fig. 5(d)). This instrument consisted 
of a pyrex cylinder with a film of nickel evaporated 
over part of the front face (fig. 8). The nickel 
films measured about 100-ohins resistance between 
contacts; the portion of the film responsible for 
this resistance was a thin strip of nickel about 
1 cm long, 1 mm wide, and 900 A thick. The 
film was connected as the active element in a d.c. 
bridge circuit. The resistance of the film was 
measured as a function of time by recording the 
output from the bridge circuit with oscilloscope- 
camera equipment. The resistance of the film 
is found to be a linear function of its temperature, 
which is taken to be the temperature of the pyrex 
surface (see appendix B). A typical temperature 
record indicated by this instrument is illustrated 
in sketch (b). Some of the actual records are 



shown in figure 9 for a variety of shock wave 
speeds from Mach number TO to 7.7. It can be 1 
seen that, to a first approximation, t tie tempera- 
ture of the glass surface jumps to a constant value 
after the reflection of the normal shock wave. 
Some irregularities appear in the output at higher 
shock speeds, apparently due to pickup of elec- 
trical charge from the hot air. The effect has 
been limited in this case by coating the nickel film 
with a thin layer of silicon monoxide. A further 
discussion of this effect and a description of the 
preparation and calibration of the films is given 
in appendix B. The temperature measurements 
are reproducible to the order of ±20 percent. 

Where the surface temperature jumps to a 
constant value, To, the heat flux into the wall can 
he calculated if the coefficients of thermal con- 
ductivity fe s and diffusivity a s are both known 



Pressure, p, atm Pressure, P, atm 


lf» 


TECHNICAL REPORT R-27- -NATIONAL AERONAUT!! S AND SPACE ADMINISTRATION 



£ 



Time, /, ^sec 


(b) 




Time, t, /i sec 
(e) 

(a) A/\ = :i.04, pi~ 1.00 atm, station 4. (b) M\ = 4.05, = 0.250 atm, station 4. 

(c) a/", = 5.44, pi ~0.002 atm, station 5. (d) = 6.72, p { = 0.010 atm, station 0. 

(e) J\/| — 7.79, pi — 0.004 atm, i tation 0. 


Figure (i. — Shock-tube pressure records. 
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(a) Pressure of air following the normal shock wave. 

(b) Pressure of air after reflection of the normal shock wave. 

Fiouke 7. — 'Measured pressures in the shock tube as a 
function of calculated equilibrium pressures. 

for tile wall material (see eq. (25)). In the 
present experiments, the following values have 
been used for the thermal coefficients of pvrex 

a s — 0 . 00 5 7 om 2 /sec ( 2 (>a ) 

k s = 0.0024 cal/cm deg see (20b) 

The^diffusivity a 0 and the integral of thermal 
conductivity <po for the air at the wall art* also a 
function of the measured surface temperature 
T q . At moderate temperatures, such as obtained 


in the present experiments, these art* approximately 
given by 

v?o = 'h2X10~ 5 7o 3/2 cal/cm sec (27a) 

a 0 = 17.7 (<po/p) cm 2 /sec (271)) 

where 7 7 0 is in degrees Kelvin anti p in atmospheres. 

Tlie quantities given by equations (20) and (27) 
permit the evaluation of the slope (dz/dy ) o from 
the measured interface temperature r ]\ (eq. (25)). 
Then the value of and thus may he ob- 
tained from the solutions shown gra])hically in 
figure 2. Since the value of or T Z} is known 

from the measured shock speed, the objective of 
the experiment is completed; that is, p is evaluated 
as a function of 7 T . 

DISCUSSION OF RESULTS 

The results from the experimental evaluation 
of ip as a function of T are shown in figure 10. 
The dependent variable is the ratio of <p to the 
value (p* for idealized air with constant specific 
heat. The latter is given by 

v?*=3.2X10 -6 T 3/2 cal/cm sec (28) 

where T is in degrees Kelvin. The ratio <pf<p* is 
approximately unity as long as the air remains in 
the molecular state. The oxygen molecules begin 
to dissociate at about 2,500° K and <p becomes 
larger than its ideal value as a result of the in- 
crease in thermal conductivity caused by the 
dissociation reaction. As the temperature in- 
creases to the point where oxygen dissociation is 
nearly complete, <p becomes about 2 l / 2 times 
larger than its ideal gas value. The theoretical 
curve shown in figure 10 was computed from 
graphical integration of the coefficients of thermal 
conductivity presented in reference 12. The 
experimental values lie somewhat above the 
theoretical curve, but the difference is hardly 
more than the* deviation in the experimental data. 
It should be noted that the function shown is for 
the particular progression of pressures obtained 
in the shock tube. Similar curves would be 
obtained at other pressure levels. For example, 
the theoretical values for <p are shown in figure 1 1 
as functions of temperature up to 15,000° K and 
for constant pressures of 0.0001, 0.001, 0.01, 
0.1, 1.0, 10 and 100 atmospheres. As temperature 
increases, the three distinct jumps in the value of 
p are due to the dissociation of oxygen, the dis- 
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Fine re 8. — Surface temperature detector. 


social ion of nitrogen, and the single ionization of 
both oxygen and nitrogen atoms. Each of these 
reactions causes a large increase in the coefficient 
of thermal conductivity (refs, (i and 12) and a 
corresponding increase in the integral, <p. 

The experimental data suggest that the theo- 
retical estimates of thermal conductivity for air, 
in which oxygen dissociation occurs, may be low 
by about 20 percent. 'Because of uncertainties in 
estimating collision cross sections (see ref. 12), 
a discrepancy of this magnitude in the theory is 
easily conceivable. The difference could also be 
caused by the presence of free electrons in excess 
of the concentrations predicted for pure air in 
equilibrium (ref. 17). According to reference 12, 
electron mol fractions {lie order of 1()~ 3 would be 
sufficient to increase the calculated thermal con- 
ductivity to the level of the present, experimental 
data. Such electron concentrations might result 
from small amounts of impurities in the shock- 


heated air. However, the inconsistencies in 
surface temperature measurements are at present 
responsible for deviations of the same magnitude 
as the lifference between the experiment and the 
thcoret cal estimates. Therefore, more precise 
experimental data will be needed before it should 
be inferred that a significant difference does exist. 

CONCLUDING REMARKS 

In calculating the flow of heat through gases, 
it is important to consider the variation in the 
coefficients of thermal conductivity and difTusivity. 
The natural parameter which replaces temperature 
as the dependent variable in such calculations is 
the integral of thermal conductivity with temper- 
ature. With this choice of the dependent variable, 
the numerical integration of the one-dimensional 
heat-ccnduction equation is straightforward if the 
initial and boundary conditions are constants and 
if the coefficient of difTusivity is assumed to be a 
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Air with equilibrium j 

oxygen dissociation ref. (12) l Theory 



0 12 3 4 5 6xl0 3 

Temperature, T , °K 

FmntE 1 0. Integral of thermal conductivity for air as a 
function of temperature; ^ (2.8) ko 7 T 0 ( Tf Tq ) 3 12 . 

known function of the conductivity integral alone. 
For ideal gases, diffusivity is approximately 
proportional to the conductivity integral. Air 
can be treated as an ideal gas up to the point 
where oxygen dissociation becomes effective. 
Over the range of temperature where the oxygen 
dissociation occurs, it is reasonable to assume that 
diffusivity is independent of the conductivity 
integral. The solutions obtained in this way 


break down at temperatures where oxygen dis- 
sociaticn becomes complete. 

The shock tube is an instrument wliich can 
produce a reservoir of high-temperature gas with 
approximately constant initial conditions. These 
initial conditions may be calculated from measure- 
ments of the shock velocity. Temperature meas- 
urements at the surface of a wall which reflects 
the shock wave can be used in conjunction with 
the solutions for heat conduction in the gas to evalu- 
ate the conductivity integral as a function of tem- 
perature. At the present stage of development, the 
experiments exhibit considerable variation, mainly 
because of uncertainty in the measurements of 
surface temperature. However, the experiments 
show the same functional behavior for the co- 
efficient of thermal conductivity as predicted by 
theoretical estimates. The agreement between 
the experiment and theory is within the accuracy 
of the measurements and within the limits of 
uncertainty in theoretical calculations of transport 
propcrfe.es. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 27 , 1959 
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Figure 11. — Theoretical estimates for the integral of thermal conductivity for air; <£* (2/Il)/; u T , (j( T; 7\,) 3 ;2 . 



APPENDIX A 

CHEMICAL RELAXATION IN SHOCK-IIEATED AIR 


It has boon assumed in this report that the 
shock-heated air reaches instantaneous equi- 
librium. It is well known that actually a finite 
number of molecular or atomic collisions are 
required to produce the equilibrium partition of 
energy (ref. 18). Thus, the gas takes a finite time 
to relax to the steady-state condition. It is the 
purpose of this appendix to discuss the question 
whether the relaxation phenomena might affect the 
heat flux experiments in the shock tube. 

At the temperatures and pressures encountered 
in the present experiments (see table IT), the 
following modes of energy can be excited in air: 
the kinetic motion of molecules and atoms, the 
rotation of molecules, the vibration of molecules, 
and the dissociation of oxygen molecules into 
atoms. No more than a few collisions are re- 
quired to produce the equilibrium distribution of 
kinetic and rotational energies. Thus, for prac- 
tical purposes, the relaxation times for these 
modes of energy are negligible. Of concern, then, 
are the vibration of oxygen and nitrogen molecules 
and the dissociation of oxygen. 

The vibrational relaxation times have been 
calculated for the present experiments using the 
measurements reported by Blackmail (ref. 19). 
The relaxation times for oxygen molecule vibra- 
tions were very short, the order of a fraction of a 
microsecond. This interval is less than the time 
resolution of the instruments involved. Nitrogen 
molecule vibrational relaxation times were some- 
what longer, the order of several microseconds, 
and these are listed in table IT. Although the 
instruments might resolve these time increments, 
tlie relaxation times are still very short compared 
to tht‘ total test interval (the order of a hundred 
microseconds or more). In addition, the energy 


transfer involved in the excitation of nitrogen 
vibrations is at most 18 percent of the total 
enthalpy. Therefore, the assumption of equilib- 
rium appears to be justified in this case. 

In the case of oxygen dissociation, the energy 
transfer by the reaction is comparatively large, 
and the corresponding relaxation times appear 
to be somewhat longer than in the ease of vibra- 
tional excitation. The oxygen dissociation re- 
laxation times have been estimated from the 
experimental results of Oaniac, Camm, Keek, 
and Petty (ref. 20). These estimates range from 
10 to 20 microseconds as shown in the last column 
of table II. These times are still somewhat 
shorter than the test, interval, and the wall tem- 
peratures appear to approach a steady value. 
In addition, recall that the measured values for 
]) 3 approach the calculated equilibrium values 
(fig. 7 1 1)) ) . Thus the assumption of equilibrium 
may be a reasonably good first approximation. 

A precise analysis of the shock reflection experi- 
ment would need to include consideration of dis- 
sociation relaxation. The gas near the wall 

would be at a higher initial temperature than the 
gas away from the wall, which would have time to 
approach the lower equilibrium temperatures 
before being influenced by the wall conditions. 
After lie reflection of the shock wave, the gas 
near t ic wall will approach equilibrium at the 
same time it is cooled by heat conduction to the 
wall. The relaxation process is equivalent to 
having a heat sink in the gas. Thus the exact 
solution requires ti e calculation of heat conduc- 
tion 1 1 rough a medium with variable coefficients 
and nonconstant initial conditions, and in the 
presence of a heat sink function distributed in 
time and space. 



APPENDIX B 


PREPARATION AND CALIBRATION OF NICKEL-FILM 
SURFACE TEMPERATURE DETECTORS 


The instrument used to measure the tempera- 
ture at the surface of the shock reflecting wall 
consisted of a thin strip of nickel evaporated 
onto a pyrex glass backing (fig. 8). Glass was 
chosen for the backing material because its sur- 
face temperature, resulting from heat conduction 
through the gas, would be higher than for most 
other solids. This is due to the fact that glass 
has a small value for the product kCp, and the 
interface temperature varies inversely as the 
square root of this product (see ref. 13). Physi- 
cally, it can be seen that a material which conducts 
the heat away from the surface slowly and which 
has a small thermal capacity per unit volume will 
come to a comparatively high surface temperature. 
Nickel was used for the film resistance element 
because it is relatively easy to evaporate with 
high purity and it has a large coefficient of re- 
sistance change with temperature. 

The nickel film is thin enough so that it is in 
equilibrium with the surface temperature of the 
glass. For example, if a film of thickness Ax is 
backed by a material having a lower value of 
k( p, the time r for the average 4 temperature of the 
film to come within 0.01 of a step input in tem- 
perature is limited by 

(A x ) 2 

T< KirtXJO 4 

where a is the diffusivity of the nickel (about 0.10 
cm 2 /sec). Thus, for a film 900 A thick, the time r 
is loss than 0.32 microsecond. This time is so 
short that the film essentially responds instantly 
to changes in the surface temperature. 

The resistances of the films were calibrated as a 
function of temperature between 20° and 90° ( 1 by 
a static method. The film, the glass slab, and t lie- 
mounting wore placed in a container immersed in 
a water bath. Thermostats and heaters in the 
bath controlled the temperature in the container 
wit inn 0.5° C. A thermocouple was imbedded in 
the metal mounting and another was placed in the 


air of the container about 2 mm directly above the 
nickel film. The calibration measurements were 
made after the output from these thermocouples 
agreed within 0.5° C. The thermal masses in- 
volved wen 4 such that about 5 minutes were re- 
quired for the film mounting and the air within 
the container to reach equilibrium with the bath. 

The resistance of the nickel films was measured 
using a standard bridge circuit and a potentiom- 
eter (in tiie shock-tube tests an oscilloscope- 
camera recorder replaced the potentiometer in 
order to achieve tin 4 time resolution required). 
It was found that if the nickel films were carefully 
prepared, the resistance was a linear function of 
temperature within the range from 20° to 90° (\ 
The deviation of the calibration data from tin 4 
least squares linear fit was about 1 percent . The 
temperature coefficient of resistivity, a , was com- 
puted from tlu 4 slope, m , of the resistance versus 
temperature data divided by the reference value 
of resistance, //o, at 20° ( 1 



or 
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It was found that the coefficient for the evaporated 
films was considerably smaller than t he coefficient 
for bulk nickel. For example, the average value 
of a for 10 stable film instruments was 0.0017 per 
deg ( ■, with the variation between films amounting 
to about 10 percent. In contrast to this, the 
coefficient for bulk nickel at 20° C is about 0.0043 
per deg (■ (ref. 21). liabinowicz, Jessey, and 
Bartseli (ref. 22) found that the same order of 
difference exists between the coefficients for 
platinum, films and hulk platinum. 

Several of the difficulties encountered in the 
preparation and use of the nickel films are note- 
worthy. First of all, the glass surface and the 
vacuum system in which it is placed must be clean 
before evaporation of the nickel is attempted. It 
was found, for example, that the nickel films could 
have a very nonlinear calibration where the 
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vacuum system had been contaminated with 
pump oil or other organic vapors. Moreover, the 
calibration was not always stable, and the values 
of a wen* generally considerably lower than 0.0017 
per deg (\ The nickel used in the evaporation 
should be quite pure also, of course. According 
to reference 22 similar precautions are essential 
in the preparation of platinum films. 

A second property of the nickel films which has 
been observed is that their calibration is not stable 
until the films have been baked at elevated tem- 
perature for an extended time. The films were 
baked at 180° C for about 1(> hours. Generally 
this resulted in a film with stable characteristics 
and one which was considerably tougher than the 
unbaked film. No attempt was made to deter- 
mine an optimum baking procedure. 

When tl ic films were used in the shock tube, 
they sometimes increased slightly in resistance 
after a test run. This might be due to erosion of 
the film by small particle bombardment. Evi- 
dence of such erosion could be seen on microscopic 
examination of the film. The resistance increase 
could also be associated with the rather large in- 
crease in film temperature after the test interval 
(see fig. 9). This increase in temperature* is prob- 
ably due to the reflection of shock waves from tin* 
helium interface and to the large increase in 
thermal conductivity and diffusivitv of the hot air 
when it becomes mixed with small amounts of 
helium. Similar increases in film resistance were 
observed in a calibration cycle after overheating. 
Euless the film was permanently damaged, the in- 
crease in resistance was less than 2 percent, how- 
ever. Films wen* calibrated both before and after 
a series of shock-tube runs, and it was found that 
the slope of resistance versus temperature was not 
affected by the change in the reference resistance. 
The value of a is a function of 11 0 , however (see 
eq. (Bl)), and it was adjusted accordingly. 

Finally, one more troublesome effect which oc- 
curred when the films were used in shock-heated 
air above 2,500° K was due to the electric charge 
picked up by the film. The effect was apparently 
the same as reported by Jahn and Weimer (ref. 23). 
The charge pickup results in spurious signals 
much larger than those due to the temperature 
when the bridge circuit is grounded. Even with 


a floating bridge circuit the spurious signals were 
the same order of magnitude as the temperature 
effect. Therefore, the nickel films used at shock 
speeds above Mach number 4.5 were coated with 
a thin evaporated layer of silicon monoxide. 
This effectively reduced the charge pickup, hut 
also increased the response time of the detector. 
A layer of SiO about 900 A thick was used over 
the nickel films in the present case. If the tem- 
perature distribution through the silicon is given 
by 

T=T a -(T 0 -T t ) erf —== 

\4a s t 

then measurements of surface temperature at the 
depth Ax from the surface will he in error by about 
Ax/\4ai. If Ax is about 10 -5 cm and the difTusivity 
a is about 0.008 cm 2 /sec for SiO, then the error at 
10 microseconds is about 5}£ percent. Tins is 
considerably less than the variation in the data 

(fig; 10- 

No attempt was made to optimize the silicon 
monoxi le coating procedure. It was found, how- 
ever, that the nickel films should he baked prior to 
application of the coating, and that the SiO should 
he* evaporated in a clean environment onto a clean 
film. Otherwise, erratic and nonlinear calibration 
of the film resulted. Even the most carefully 
coated ilms varied from one another in the value 
of a by as much as 20 percent, about twice* the 
variation observed in the uneoated films. This 
variation is accounted for in the calibration, of 
course, but it may be an indication of some addi- 
tional instability in the film. 

In conclusion, the nickel film surface tempera- 
ture de Lectors developed for the present experi- 
ments v ere found suitable for preliminary investi- 
gations of the heat flux through gases. However, 
they are at present the weakest link in the experi- 
ments, end most of the variation in the data (fig. 
10) is attributable to instabilities in this instru- 
ment. Temperature detectors of greater stability 
and dui ability, of more uniform calibration, and 
free from charge accumulation or shunting effects 
are nee< ed for more precise determination of the 
thermal conduction properties of high-temperature 
gases. 
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